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The present paper suggests the continuum theory of Tkachenko modes in a rotating 2D Bose- 
Einstein condensate taking into account density inhomogeneity and compressibility of the conden- 
sate. The problem requires solution of coupled hydrodynamic equations for vortex and liquid motion 
with proper boundary conditions, which were derived for the condensate described by the Thomas- 
Fermi approximation. Compressibility becomes essential at rapid rotation with angular velocity 
close to the trap frequency. The theory is in a reasonable agreement with experimental observation 
of Tkachenko modes. 
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Tkachenko mode 1] is a transverse sound wave in a 
vortex lattice, which exists in the ground state of a ro- 
tating supcrfluid. It has attracted a lot of attention in 
vortex dynamics of laboratory and astrophysical super- 
fluids. Observation of this extremely soft mode (of order 
or less than a few Hz) in laboratory superfluids was very 
difficult because even slight pinning of vortex ends shades 
the contribution of vortex shear rigidity transforming 
the Tkachenko wave into a classical inertial wave [J. 
However, discovery of Bose-Einstein condensate (BEC) 
and possibility to rapidly rotate it made observation of 
Tkachenko modes more feasible, which resulted in clear 
experimental detection of them |2j| . As well as the Kelvin 
mode, also possible in the BEC [5| , the Tkachenko mode 
is a direct manifestation of the quantum vorticity since 
it depends on the circulation quantum. 

In a number of aspects external conditions for 
Tkachenko waves in BEC are essentially different from 
those in "old" superfluids. First, in a rapidly rotating 
BEC one cannot consider a superfluid to be incompress- 
ible. Effect of finite compressibility was investigated the- 
oretically in the past [2| since the effect is crucial for 
the long-wavelength limit of the Tkachenko mode. But 
at that time it was considered as a theoretical curiosity 
since in order to reveal it one needed containers about 
a few hundred meters diameter. In the BEC case the 
situation is essentially different because in contrast to a 
strongly interacting bose liquid as He II, BEC is a weakly 
interacting bose gas with very low sound speed and very 
high compressibility. Importance of high liquid com- 
pressibility for Tkachenko waves in BEC was pointed out 
by Baym [J], who rederived the spectrum of Tkachenko 
waves in a compressible liquid [|| known from Ref. and 
compared it with the experiment 0] ■ Another important 
feature of rotating BEC, also connected with its high 
compressibility, is that the liquid density is essentially in- 
homogeneous. This feature was taken into account in the 
theory by Anglin and Crescimanno within the frame 
of the continuum theory, which replaces a discrete vortex 
lattice by a continuous medium, like the elasticity the- 
ory for atomic crystals. But they neglected liquid com- 



pressibility, while a proper comparison with the experi- 
ment requires a theory taking into account the both fea- 
tures, compressibility and inhomogeneity. Recently the 
Tkachenko mode in rotating BEC was investigated nu- 
merically with solving the equations of Gross-Pitaevskii 
theory (mean-field theory) 0, @ ■ The numerical results 
well agreed with the experiment. But numerical calcu- 
lations, whatever powerful and useful they can be, do 
not exclude necessity to develop the analytical approach, 
as far as it possible, since the latter gives analytical and 
scaling properties and therefore provides a deeper insight 
into physics of the phenomenon. In the experiment there 
are good conditions for application of a continuum the- 
ory since usually the relevant length scale (the condensate 
size) essentially exceeded the intervortex distance. The 
goal of the present paper is to suggest the continuum 
theory of Tkachenko modes in a rapidly rotating BEC 
taking into account liquid compressibility and inhomo- 
geneity. The main challenge for the theory was to formu- 
late proper boundary conditions for oscillating BEC. 

Let us remind the equations of motion of a homoge- 
neous compressible superfluid at T = (no normal com- 
ponent) in the rotating coordinate frame Q: 



/<„ ( — + 2n x v L 



2v7 / 



p K X (vl - V) = 



u) - Ai 



(1) 



(2) 



(3) 



Here p'(f) is the oscillating component of the liquid 
mass density p(r) = p + p'(r) around the equilibrium 
homogeneous density p , u is the vortex displacement, 
Vl = duj dt is the vortex velocity, v is the liquid velocity 
averaged over the vortex-lattice cell, fi is the angular ve- 
locity, c s is the sound velocity, k — h/m is the circulation 
quantum, and ct = y kH/8it is the Tkachenko-wave ve- 
locity. Assuming that ct <C c s and crk <C f2 one receives 
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for the spectrum of plane waves oc exp(ifc • r 
dispersion relation: 
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which yields the gapped sound mode 
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and the soft quantum mode (see Ref. [l(J and Eq. (4.73) 
in Ref. @) 
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So for the soft (Tkachenko) mode compressibility is essen- 
tial in the long- wavelength limit k <C Q/c s . It transforms 
the Tkachenko wave with the sound spectrum lu = c^k 
to a softer mode ui oc k 2 . In the Tkachenko mode motion 
of both the vortex lattice and the liquid has an elliptic 
polarization but longitudinal components parallel to the 
wave vector k are very small. The transverse components 
(normal to k) of the vortex lattice and the liquid are close 
one to another: VLt ~ Vt- 

Generalization of these equations onto an inhomoge- 
neous liquid with equilibrium density po (r ) , which varies 
in the plane normal to the rotation axis, is straightfor- 
ward. The goal of this paper is to consider only axisym- 
metric Tkachenko eigenmodes in a axisymmetric rotating 
BEC. Thus we write the equations of motion in the po- 
lar system of coordinates for the monochromatic mode 
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We excluded all variables except for the tangential (az- 
imuthal in polar coordinates) component of the velocity 
Vt ~ and the radial component of the liquid velocity 
v r . The latter, though being much smaller than vt, is 
crucial for the compressibility effect. For a weakly in- 
teracting bose gas c 2 is proportional to the density po. 
Therefore the ratio c 2 /p Q can be replaced by its value 
cg(0)/p o (0) in the BEC center r = 0. 

Let us remind the results for Tkachenko eigenmodes in 
an incompressible (c s — > oo) homogeneous (po = const) 
liquid in a cylindric container. Then the radial liquid 
velocity v r vanishes, and Eq. Q is the Bessel equation 
with a solution v t oc J\(kr). The eigenvalues of the wave 
number k are determined from the boundary condition 
at r = R, where R is the radius of the vortex-lattice 
sample, which normally is close to the container radius. 
In the ideal axisymmetric case without any interaction 
of vortices with lateral walls the radial flux of the trans- 
verse momentum of the vortex lattice, which is given by 



the corresponding component of the lattice stress tensor, 
must vanish. This gives the boundary condition 
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The flux of the angular moment into the liquid is propor- 
tional to the same stress-tensor component and there- 
fore vanishes also. Equation © leads to the condition 
JzQcR) — imposed on the wave numbers k. Tradition- 
ally in the papers on BEC they scale the Tkachenko- mode 

frequencies by the frequency fib/ 1 R, where b — \/ k/*</3Q 
is the intervortex distance 
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Then the first two roots kR — 5.14 and 8.42 of the equa- 
tion J2{kR) = give the first two reduced eigenfrequen- 

cies of the Tkachenko mode to = J \/3/87rfci? = 0.263A:i? 
with the ratio 8.42/5.14=1.64. Another interesting case 
is the limit of strong interaction of vortices with rough 
lateral walls, when v t (R) = 0. Then the eigenfrequencies 
are given by two first roots of the equation Ji(kR) = 
kR = 3.83 and 7.02 with the ratio of the two lowest 
eigenfrequencies 7.02/3.83=1.83. 

Now let us consider an inhomogeneous (but still incom- 
pressible) liquid. The analysis will be applied to the ro- 
tating 2D BEC in pancake geometry with the density pro- 
file po = p(0)(l — r 2 /R 2 ) determined from the Thomas- 
Fermi approximation [ll| . Here R = \/2c s (0) / yJio\ — O 2 
is the BEC radius (Thomas-Fermi radius), c s (0) is the 
sound velocity at the symmetry axis r — 0, and u>± is 
the trap frequency, which characterizes the curvature of 
the trapping parabolic potential Since the momen- 
tum flux is proportional to po and the latter vanishes at 
r = R, it looks that the flux through the boundary van- 
ishes independently on whether the boundary condition 
Eq. © is satisfied or not and the latter is unnecessary. 
But this is not the case as one can see from solution 
of the equation of motion close to r = R by expansion 
in small (R — r)/R. With the relative error of the or- 
der of (R — r) 2 /R 2 the general solution of Eq. JJJ) is 
Vt ps r[C\ + C*2 ln(i? — r)], where C\ and C2 are arbitrary 
constants. The component oc C'2 gives a divergent contri- 
bution to the displacements and a finite contribution to 
the stress tensor despite the factor p oc R— r. This would 
violate the conservation law for the angular moment and 
this component should be absent. This requirement is 
satisfied if the boundary condition Eq. © takes place. 
Therefore the condition Eq. is relevant also for an 
inhomogeneous liquid with po(r) — > at the border. The 
equation of motion for an inhomogeneous liquid has no 
analytic solution in the whole interval of r but it is not 
difficult using Mathematica to find the eigenmodes. First 
two of them are u> — 1.428 and 2.327 with the ratio of 
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the two lowest eigenfrequencies 2.327/1.428=1.63. This 
agrees with the results of Ref. 0- 

According to Eq. JBJ the compressibility effect be- 
comes important if k is of the order or less than Sl/c s . 
Since the eigenvalues of k are of the order of 1/R this 
yields the condition for a strong compressibility effect: 

nR/ Cs (o) 
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SI 2 ^> 1. Thus at rapid rotation 



of the BEC with angular velocity SI close to the trap 
velocity ui±. liquid compressibility should be taken into 
account. If we want to make the quantitative analysis for 
an incompressible liquid, we should solve the system of 
two coupled second-order differential equations Q and 
© with proper boundary conditions. One of them re- 
mains to be Eq. |JS}, but we need also the boundary 
condition imposed on the radial liquid velocity v r . We 
use the arguments similar to those used for derivation of 
Eq. The total mass balance requires that the ra- 

dial mass current po{r)v r (r) at the BEC border r = R 
vanishes. Indeed, even though the BEC border is mobile 
because of the oscillation, the total mass, which can be 
transferred through the equilibrium border r = R, is a 
second-order quantity with respect to an oscillation am- 
plitude. The radial mass current vanishes for any finite 
radial velocity v r (R) since p(r) — > at r — > R. But some 
condition is required to provide v r (R) to be finite. In 
order to derive this condition we solve Eq. JHJ at r w R 
by series expansion [neglecting terms ~ (R — r) 2 ]: 
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The divergent component oc C2 gives a finite mass flow at 
the border and should be ruled out. Taking a derivative 
from this expression and excluding the constant C\ we 
receive the boundary condition imposed on v r : 
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The left-hand side of the boundary condition (|12f> is the 
divergence V-u of the liquid velocity in polar coordinates. 
The divergence remains constant at the boundary of the 
BEC sample, but the density po(R) vanishes there. As a 
result, the mass continuity equation yields that 
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Thus at the sample boundary the whole variation of the 
density originates from the oscillation of the boundary 
with the velocity v r (R). The same boundary condition is 
valid for other compressional modes of the BEC [I2I Il3] | . 

For better understanding of parametric dependences of 
the eigenfrequencies it is useful to introduce the dimen- 
sionless radius z = r/R and to scale the radial velocity, 
introducing the velocity u: v r = iucT/c s (0). Then the 
system of equations is: 

sq(l-z 2 )u(z) = q 2 (l- z 2 )v t (z) 



l_d_ 

z 2 dz 



(1 - z 2 )z 3 



dz 



v t (z) 



d f 1 d\(l - z 2 )zu(z)] 

sqv t (z) = — i 

dz [ z dz 



(14) 



(15) 



where q determines the reduced frequency, u> — 0.263g , 
and the parameter 
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characterizes the effect of compressibility. In the dimen- 
sionless form the boundary condition (|12fl is du(l)/dz + 
u(l) = sqv t (l)/4. 

Let us discuss now restrictions on validness of the pre- 
sented analysis. In order to observe the strong compress- 
ibility effect s ^> 1 one should rotate the BEC rapidly, 
at SI ~ But if SI approaches to uj± close enough at 
fixed number of atoms the condition of small ratio ct / c s 
is violated. Since the vortex core size is £ ~ n/c s , the 
large ratio ct/c s ~ £/b means that the vortex cores start 
to overlap. This signals that we approach to the critical 
angular velocity, similar to the upper critical magnetic 
field H C 2, when rotation should suppress the superfluid 
order parameter. When £ ~ b the properties of the BEC 
become essentially different from those assumed in the 
present work (see Ref. Q and references therein). How- 
ever with a number of atoms large enough one can reach 
very high values of s without violation of the condition 
<C 1 UM- Another assumption of our analysis was 
solid body rotation of the BEC with the constant vortex 
density given by the Feynman formula n v — 2Sl/n. In a 
inhomogeneous liquid this formula is not exact. Correc- 
tions to this formula due to forces from density gradients 
were found in Refs. Q andE But these corrections 
are essential only at the distance of the order b from the 
sample boundary 0] and therefore are not essential for 
determination of eigenfrequencies until R~^> b. 

In order to find the eigenfrequencies we determined 
numerically (at fixed s) eigenvalues of q at which Eqs. 
ifHj) and fT5|) have a solution with proper boundary con- 
ditions. In addition to the boundary conditions at the 
BEC boundary discussed above, there are standard con- 
ditions in the BEC center that the velocities are not 
divergent, namely, v r oc v t oc r. If compressibility is 
important any reduced eigenfrequency is not a number 
but a function of s: ibi = fi(s). At s — > the func- 
tions fi(s) are constant, but at large s the functions fi(s) 
are inversely proportional to s. The numerical solution 
yields that at large s the two first eigenfrequencies are 
u>i = 7.17/ s and G)% = 16.9/s. Figure ^ shows the nu- 
merically found first eigenfrequency u>i plotted as a func- 
tion of n/^Lol - SI 2 = s/2V2 (solid line). The black 
squares show experimental points plotted in dimen- 
sionless variables by I. Coddington. They were obtained 



4 




1 2 3 4 5 6 7 

ft 

^ j - o, 2 

FIG. 1: Comparison between the theory (solid line) and the 
experiment (black squares). 



for various parameters, but collapse on the same curve, as 
expected from the present analysis. Quantitative agree- 
ment between the theory and the experiment looks quite 
good. It becomes worse at larger fi/ — Q 2 , since the 
growth of the parameter £/& in this area makes the theory 
less accurate (see above). Coddington al. Q measured 
also the ratio of the two first frequencies u>2/u>i = 1.8 
at il/uj± = 0.95, which corresponds to s = 8.61, The 
present theory predicts the ratio ^2/^1 = 2.09. 

In summary the continuum theory of Tkachenko modes 
was developed, which takes into account the liquid in- 
homogeneity and finite compressibility. Application of 
the theory to a rapidly rotating Bose-Einstein conden- 
sate required formulation of proper boundary conditions 
at the condensate border where atom density vanishes. 
The theory is in good agreement with observations of 
Tkachenko modes in the rapidly rotating Bose-Einstein 



condensate. 

I thank Ian Coddington who prepared the experimen- 
tal part of Fig. ^ I also appreciate very much an interest- 
ing comment by Marco Cozzini, Lev Pitaevskii, and San- 
dro Stringari concerning comparison of the present anal- 
ysis with that based on the compressibility sum rule Q. 
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